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TVj = 1 QCD simulation with improved gauge action at finite temperature 
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We study the deconfinement transition of Nf — 1 QCD by the hybrid Monte Carlo algorithm with Wilson 
fermions. We calculate the Polyakov loop, its susceptibility and Binder cumulant and use the \ 2 method to locate 
the phase transition point. Our results are similar to the previous results obtained by the multiboson algorithm. 
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1. Introduction 

Liischer proposed the multiboson algorithm^ 
2] in which a polynomial approximation to the 
inverse of the fermion matrix is used. This 
algorithm is considered to be an alternative 
one to simulate dynamical QCD. Borigi and de 
Forcrand 3J noticed that the polynomial approx- 
imation can be used to simulate odd flavour 
QCD. Then the multiboson algorithm is used 
to study the deconfinement transition in one 
flavour QCD |3]. The original hybrid Monte Carlo 
(HMC) algorithmic is limited for even flavour 
case. However the polynomial approximation can 
be also applied for the HMC 13 , which means that 
one can also simulate odd flavour QCD by the 
HMC0. We investigate the phase structure of 
one flavour QCD by the HMC and check whether 
different algorithms give the same results. We 
calculate Polyakov loop, its susceptibility and 
Binder cumulant and use the x 2 method |S] to lo- 
cate the precise position of the phase transition. 
Mainly we present results from the Wilson gauge 
action, and also give preliminary results from the 
DBW2 gauge action]^. 



2. Nf = 1 algorithm 

The partition function of Nf = 1 QCD is given 
by Z = J [dU]detDexp(-S g [U]). One can ap- 
proximate the inverse of I? by a polynomial^ as 



1/D w P m {D) = TT(-D — Zk), where Zk are the 



fc=i 

roots of the polynomial P m (D). For the even de- 
gree m, det-D can be rewritten as 

detD = C n det(T*(D)T n (D))-\ (1) 

where T n (D) = Ut=i( D - Z M -i), and C n = 
det(DT+(D)T„(L»))" 1 is a correction factor to the 
polynomial approximation. det(T^(D)T n (D)) 
can be expressed as, 

detiTl^TniD))- 1 

\d^][dcj>]eM-^Ti{D)T n {D)cj>). 



(2) 



Then one can define Nf = 1 Hamiltonian for 
HMC algorithm as, 

i 

where Pi are momenta which are conjugate to link 
variables. ^Tj l (D)T n (D)(f> is manifestly positive, 
which is essential to the implementation of the 
HMC. The correction factor C n is used to make 
the algorithm exact [7]. 



a. X 

In this study we use x 2_m ethodjBJ to find crit- 
ical j3 for each k. By differentiating the singular 
part of the free energy density on 7V S 3 x N t lattice 
and expanding these equations with respect to x, 
we obtain the leading N s behavior at x = 0, 



ln(L)(x = 0,N s ) = c&lnJV.+ci, 



(4) 
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0.10 
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0.12 


24/0.86 15K 


32/0.75 20K 


40/0.62 14K 


0.14 


32/0.86 10K 


40/0.74 14K 


50/0.61 12K 



Table 1 

n and Acc are the number of degrees of polynomial and acceptance ratio of the HMC respectively. The 
step size dr is set to 0.05 for all simulations. To obtain the best performance of the HMC one should 
tune dr so that the acceptance is 60 — 70%^D]. In this study we have not done such tuning. 



lnx(z = 0, N s ) = colnJV a + ci, 
B 4 (x = 0,N s ) = co + Cl N-». 



(5) 
(6) 



where (L), Xi and x are Polyakov loop, its 
susceptibility, Binder cumulantjllj and reduced 
temperature (T — T c ) /T c respectively: 



pected to behave as V a . The behavior is similar 
to the previous results^], i.e. a decreases as k 
increases. However the precise values of a are 
slightly different from the previous ones, which 
may indicate that our statistics is not enough to 
conclude. 



(7) 



X = N^({L 2 ) - (I/) 2 ), Xv = N!(L 2 ), (8) 
(L A ) 



B± = 



{L 2 Y 



- 3 



(9) 



The critical point can be obtained at the point 
where a fit to the leading A^-behavior has the 
least x 2 0- 

4. Simulations 

We perform numerical simulation using the 
hybrid Monte Carlo (HMC) with polynomial 
approximation 7 . We use the Wilson and 
DBW2 9 gauge actions. Lattice gauge action is 
given by 

S g [U] = ^((l-8 Cl ) ^ ReWn+d ^ ReW 12 )(10) 

where W\\ and W\i stand for 1 x 1 and 1x2 
Wilson loop. For Wilson action the coefficient is 
set to ci = and a = -1.4089 for DBW2 action. 

To obtain values of observables at different /3 
from a single simulation, we use the re-weighting 
method^]- Statistical errors are evaluated by 
the Jackknifc method. 

5. Results 

In Fig^ we show the volume dependence of 
peak of Polyakov loop susceptibility, which is ex- 
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Figure 1. Volume dependence of peak of Polyakov 
loop susceptibility. 

We summarize results for (3 C from \L\, Xv and 
B4 in Tablc[21 We estimated the critical (3 using 
the x 2 - me thod. The data were fitted to Eq.lQJ- 
Eq. (JBJ ■ For Eq. (JBJ u was fixed to 1 in accordance 
with 8 . As an example in FigEl x 2 /Nfu curve 
at k — 0.12 is shown. Here Nf it is the degrees of 
freedom for the fitting. We determine f3 c at the 
minimum point of the x 2 /Nfu curve. 

/3c L ', /?c" an d Pc 4, should coincide with each 
other. However we observed that at k — 0.05, 0.10 
and 0.14 (3^ is slightly different from others. 
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Table 2 

Preliminary results of critical (5 from various ob- 
servables. \ P* v anc ^ Af 4 stand for critical (3 
from Polyakov loop, \ v and B4 respectively. 
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Figure 2. x 2 /N f lt curves at re = 0.12(/3 = 5.630). 
Nfu = 2. 

This might be due to less statistics. To deter- 
mine P c precisely more statistics is needed. 

The results with DBW2 gauge action ( 
quenched and n = 0.1 ) are shown in FigOl We 
are currently performing simulations at other re. 
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